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Abstract. We present an analysis of ASCA X-ray data of 
SNR RCW 86. There appears to be a remarkable spectral 
variation over the remnant, indicating temperatures vary- 
ing from 0.8 keV to > 3 keV. We have fitted these spectra 
with non-equilibrium ionization models and found that all 
regions are best fitted by emission from a hot plasma un- 
derabundant in metals (< 0.25 solar), but in some cases 
fluorescent emission indicates overabundances of Ar and 
Fe. The ionization stage of the metals appears to be far 
from equilibrium, at some spots as low as log(net) (m~^s) 
~ 15.3. We discuss the physical reality of the abundances 
and suggest an electron distribution with a supra-thermal 
tail to alleviate the strong depletion factors observed. We 
argue that RCW 86 is the result of a cavity explosion. 

Key words: Atomic Processes; ISM: abundances - ISM: 
individual objects: RCW 86 - ISM: supernova remnants - 
X-rays: spectroscopy 



1. Introduction 

RCW 86 (also G315.4-2.3 and MSH 14-63) was for some 
time regarded as the historical remnant of the supernova 
(SN) of AD 185 (Clark & Stephenson 1977). Recently, 
however. Chin & Huang (1994) presented convincing evi- 
dence that the event recorded by the Chinese was not a SN 
at all. This makes RCW 86 just one of many other galac- 
tic supernova remnants (SNRs). Nevertheless, the ASCA 
X-ray spectra show that RCW 86 has some very special 
properties. 

The radius of the remnant is ^^22'. Distance estimates 
were usually based on the notion that RCW 86 was SN185. 
The most recent (kinematic) distance estimate is 2.8 kpc 
(Rosado et al. 1996). 
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2. The data 

Our analysis is mainly based on archival ASCA data 
(Tanaka et al. 1994). RCW 86 was observed by ASCA 
in the PV phase on August 17, 1993. Three pointings 
were used. Two of the pointings (together 14 ks) cover 
the bright boomerang shaped area in the Southwest of 
the remnant. The third observation (^ 14 ks) covers the 
Northeastern part of the remnant. The Northwestern part 
of the remnant is not covered. The area covered by ASCA 
is indicated in the ROSAT PSPC image Fig. 0. 

For our analysis we used spectra from all four instru- 
ments on board ASCA (i.e. SIS O&l, CIS 2&3). The spec- 
tral resolution of the SIS instruments is better than that 
of the CIS (2% vs. 8% FWHM at 6 keV), but the latter 
have the advantage of better sensitivity at high energies. 
SIS observations were made in 4 CCD faint and bright 
modes. 

We reduced the data with the NASA/HEASARC ftools 
package v3.6 using standard screening criteria. Back- 
ground spectra were extracted from the standard back- 
ground data available at the HEASARC archive. The re- 
sponse matrices for the SIS spectra were generated with 
sisrmg, adding matrices of the various chips with the ap- 
propriate statistical weights and auxiliary response files. 

The angular resolution of the ASCA instruments en- 
sures that potential contamination by the galactic ridge is 
negligible (Claas et al. 1989, Kaastra et al. 1992). 

We also analyzed archival ROSAT PSPC data of 
RCW 86. Three pointings focus on the North and one on 
the South of the remnant. The PSPC has the advantage 
of having a better spatial resolution (FWHM ^ 30"), than 
ASCA (FWHM ~ 60"), but the spectral resolution is poor. 
PSPC data were reduced with the ftools package as well. 

3. Spectral analysis and interpretation 

Spectral modeling. As a starting point for looking for spec- 
tral variation over the remnant we made hardness ratio 
maps for all pointings and instruments. Rather than show- 
ing all ASCA ratio maps, we present here the ROSAT 
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Fig. 1. The ROSAT PSPC image of RCW 86 in the energy range 0.1-1 keV. The image is composed of three observations, 
with a total integration time of 9ks. The image was smoothed with a median filter of 7 x 7 pixels, the pixel size being 15". The 
grey scale is such that one can see a clear distinction between the faint emission (the lightest grey) and the bright shell already 
known from the Einstein (Pisarski et al. 1984) and EXOSAT missions (Claas et al. 1989). Note that the satellite pointing was 
such that the Southwestern region was observed with less sensitivity than the Northeastern region. The three squares indicate 
the areas covered by the ASCA SIS instruments. 



PSPC ratio map (Fig. Note that with the spatial res- 
olution of ASCA these regions were less distinct than in 
Fig. ^. So the spectral differences shown in Fig. ^ may in 
reality be larger. Because the ratio map does not give a 
good impression of the morphology of the remnant, we also 
present a PSPC image of RCW 86. Fig. |] shows that the 
X-ray emission extends beyond the main shell. This was 



A color version and other images are available at 
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not seen in the HRI image by Pisarski et al. 1984. In par- 
ticular the Northwest of the remnant shows faint emission 
outside the main shell. This cannot be an effect of the 
point spread function or interstellar scattering, because 
then one would expect a similar effect for the Northeast 
of the remnant, where the main shell is as bright as in the 
Northwest. Furthermore interstellar scattering takes place 
preferentially at lower energies (the photon cross section 
scales as see Predehl & Schmitt 1995), whereas Fig. 
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Fig. 2. The PSPC ratio map, based on the same observations as Fig. ^ The images are composed of three observations, with 
a total integration time of 9ks, The images were extracted in the energy ranges 0.1-1 keV and 1-2.1 keV and smoothed with a 
median filter of 7 x 7 pixels, the pixel size being 15". A higher ratio indicates a harder spectrum. The color scale is such that 
they represent our choice for selection of spectrum well. 



shows that the faint emission in the Northwest is rather 
hard. Finally note that the Ha image (Smith 1997) and the 
radio surface brightness profile (Pisarski et al. 1984) show 
also emission outside the main shell. For Ha the emission 
seems to come from breakouts from the main shell. 

The ratio map shows that we can identify two dis- 
tinct regions in the Southwest (SW). For obvious reasons 
we will call spectra from these regions SW Hard and Soft. 
The SW Soft region correlates very well with the Ha emis- 
sion (Smith 1997), whereas the SW Hard region shows a 
rough correlation with the radio map (Kesteven & Caswell 



1987). The SW Soft region is essentially the same as the 
bright boomerang-shaped shell in Fig. |l|, although the SW 
Hard region has some overlap with this shell, suggesting 
a projection effect. 

In the fainter Northeastern (NE) part of the remnant 
we also find hardness variations. The spectra of the hard- 
est region we will tag NE Hard. We extracted spectra East 
of this, including the bright arm (NE Soft) and more to the 
center of the remnant (NE Center). As Fig. || shows the 
spectra are distinctively different, but all have a charac- 
teristic peak arising from helium-like Ne around 0.9 keV. 
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Table 1. Plasma parameters for one and two temperature fits to the ASCA spectra. The abundances are number abundances 
in fractions of solar abundances (Anders & Grevese 1989). Numbers in brackets give the 90% confidence range (Ax = 2.7). 
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N 





(0-0.04) 





(0.0-0.01) 


0.02 


(0-0.2) 


0.03 


(0-0.3) 





(0-0.03) 


O 


0.09 


(0.06-0.12) 


0.014 


(0.008-0.021) 


0.01 


(0.02-0.1) 


0.012 


(0-0.06) 





(0-0.002) 


Ne 


0.25 


(0.21-0.28) 


0.078 


(0.065-0.086) 


0.08 


(0.10-0.25) 


0.08 


(0.04-0.5) 


0.04 


(0.03-0.05) 


Mg 


0.12 


(0.10-0.15) 


0.06 


(0.04-0.07) 


0.07 


(0.04-0.14) 


0.13 


(0.03-0.5) 


0.03 


(0.0-0.09) 


Si 


0.19 


(0.14-0.24) 
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(0.10-20) 
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(0.07-0.24) 
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(0-0.7) 


0.3 
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Ar 
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Fe 


0.16 


(0.14-0.20) 
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(0.09-0.16) 
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(0-> 10) 
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A^H (lO^^cm^^) 
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2.7 


(1.8-4.2) 


2.0 


(0.5-3.5) 


1.3 


(1.0-1.7) 




4.7 




1.7 




1.2 
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1.7 





The NE spectra are in general harder than the SW spec- 
tra. The correlation between faintness and hardness (the 
softest emission conies from the brightest regions) sug- 
gests that the shock has decelerated in the densest region 
compared with the fainter regions such as the Northeast. 

As noted above ASCA has a fairly large point spread 
function with broad wings making that the half power di- 
ameter is 3' (Tanaka et al., 1994, Serlemitsos et al. 1995). 
One might wonder how this affects the spectra in Fig. |[ 
Most serious is the contamination of the SW Hard spec- 
trum by emission from the SW Soft region for energies 
around 1 keV. For higher energies (> 2 keV) it is the other 
way around. As an example, based on the PSPC image and 
spectra we estimate that the contribution of the emission 
from the soft region on the spectrum of the hard region is 
less than 25% at 1 keV. Such differences are much smaller 
than the overall differences between both spectra (Fig. ||). 

We fitted the spectra from all regions with one and 
two temperature component NEI models using the spec- 
tral code SPEX (Kaastraet al. 1996). The free parameters 
were: the emission measure (EM) UeTiuV/cP (where d is 
the distance, nn and rie, the hydrogen and electron num- 
ber densities and V the volume occupied by the plasma); 
the ionization parameter Ugt (t being the time since the 
plasma was heated); the electron temperature T^; the in- 
terstellar column density N-^ and the abundances of N, 
O, Ne, Mg, Si, S, Ar and Fe. Since in some cases the 
Fe abundance can be estimated from both the Fe-L and 
Fe-K emission, we left the Fe abundance of the hot (i.e. 
Fe-K emission) and cool component free. In most cases 
the temperature of the hot component was not very well 
determined; in order to prevent temperatures much in ex- 
cess of what was found by Ginga (Kaastra et al. 1991) we 
set an upper hmit to the temperature of kTe = 5 keV. We 
have used the interstellar extinction model of Morrison & 
McCammon (1983). Table |l| hsts the results. 



Abundances. There are two remarkable things to be noted 
about the best fitting models. First, the Uct values are 
very low, in some cases as low as 2 10^^m~'^s, one of the 
lowest value ever reported {uet values a factor ten lower 
were recently reported for some regions of the Vela SNR 
based on ROSAT PSPC spectra, Bocchino et al., 1997). 
Second, the abundances are peculiar. In general, the abun- 
dances are low, except for some cases of Ar and Fe. The 
high Fe abundance is not simply an artefact of the model, 
which is clear from Fig. ^. It shows that the centroid of 
the Fe-K complex is at 6.4 keV, indicating a fluorescent 
origin. The Ar fluorescent emission at 3 keV is present in 
both the SIS (Fig. H) and the CIS spectra. Note that flu- 
orescence caused by irradiance is very unlikely in view of 
the large column densities and the large radiation field 
needed. So are the low abundances real? What could pos- 
sibly be the origin of the abundance pattern? An approach 
to these problems is to analyze the assumptions made in 
fitting the spectra. 

In the spectral code abundances are essentially deter- 
mined by the line to continuum ratios. The continuum 
is supposed to be mainly bremsstrahlung from electron- 
ion collisions, assuming a Maxwellian electron energy dis- 
tribution. Furthermore, the assumption is made that the 
plasma can be approximated by a simple one or two com- 
ponent plasma model assuming constant temperature, al- 
though in reality temperature and density gradients must 
be present. For most SNRs such a model works reason- 
ably well. A possible extra source of continuum radia- 
tion can originate from synchrotron radiation arising from 
shock accelerated electrons with energies of the order of 
a few TeV. Such a model has been proposed for SN 1006 
(Koyama et al. 1995). However, in our case this explana- 
tion does not work. First of all, the SW Soft spectrum 
is too soft to be dominated by a synchrotron component 
which is usually assumed to have a power law behaviour. 
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Fig. 3. SIS spectra of various regions together with NET model 
spectra (solid lines). From top to bottom; SW Soft (xlO), SW 
Hard, NE Soft, NE Center (xO.l), NE Hard (x0.005). 



The second reason is even more compelling: a power law 
component makes it even harder to understand why we 
see such a relatively strong Fe-K complex around 6.4 keV. 

Another way to deal with under abundances, dust de- 
pletion, is very unlikely in this case because the abundance 
pattern wc find is not what one expects it to be in case of 
dust depletion. For instance Ne is under abundant, whereas 
it cannot be depleted in grains. Nevertheless, although it 
does not explain our overall results, dust depletion may 
give an additional effect for some metals (O, Mg, Si, Fe). 

One might also think that an increased He/H ratio 
with respect to the solar value results in a stronger con- 
tinuum emission and thus in a smaller line to continuum 
ratio. However, this is not the case. To see this, consider 
two plasmas, one with a solar He/H ratio and one for 
which part of the H has been used to mak(^ He (mostly 
^He). Because 4 H ions provide more electrons than one 
He ion, the plasma with the solar He/H ratio will emit 
more bremsstrahlung continuum. 

So we probably need a modification of the standard 
model. A possibility is that the electron distribution is not 
Maxwellian. The continuum emission may partly result 




6.5 

Energy (keV) 

Fig. 4. The Fe-K complex in the South as observed by GIS2 
and GIS3 (combined in one spectrum). The centroid of the 
complex is 6.3 ± 0.1 keV consistent with a fluorescent origin. 



from a (non-thermal) tail to a relatively cool ( ^ 0.05 keV) 
electron distribution. As a result, metals should be in a 
lower ionizaticm stage than implied by the model fits. Such 
a model does at least qualitatively explain why fluorescent 
emission indicates different (i.e. high) abundances. 

We have tried to model such a distribution by calculat- 
ing the ionization stages of O for a mix of a high (~ 5 keV) 
and low temperature electron distribution (~ 0.02 keV). 
This crude model indicates that we can hide an apprecia- 
ble amount (40%) of metals as long as riet is low (e.g. 
^510^^m~^s for a cool to hot ratio of 4 to 1). Note 
that an appreciable amount of cool electrons is not un- 
likely in the light of recent UV observations of SN 1006 
(Laming et al. 1996), which indicated Te/Tio„ ^ 0.05. As 
an indication: a shock velocity for RCW 86 of 700 km/s 
(Long &: Blair 1990) without electron-ion equilibration 
gives kTe = 0.0005 keV; Laming et al.'s result suggests 
kTe ^ 0.05 keV. Clearly the effects of a non-Maxwellian 
electron distribution need further modeling. It is further- 
more unclear why all the spectra of RCW 86 share so many 
characteristics, whereas they are taken from very different 
regions. 



Up to now we did not address the issue how projec- 
tion or contamination effects may be responsible for the 
unusual ASCA RCW 86 spectra. In our opinion it is very 
hard to explain the near absence of line emission solely 
by projection effects without invoking an unconventional 
interpretation. In a model that uses projection effects the 
observer should see a superposition of a spectrum with 
a "normal" line to continuum ratio and a spectrum with 
an almost complete absence of lines. So the difficulty of 
how to suppress line emission remains unaltered in such 
a model. In fact we cannot exclude that the ASCA SIS 
spectra are superpositions of two opposite cases, but we 
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Fig. 5. A sketch of our interpretation of the structure of RCW 86 as a cavity explosion. We have drawn the blastwave here as 
a sphere, but in reality the shape of the blastwave will be affected when it hits the cavity wall. Bright are those regions where 
the blastwave has hit the cavity wall. 



want to point out that it is remarkable that such extremely 
smooth spectra exist. 

Evidence for a cavity explosion. Pisarski et al. 1984 and 
Claas et al. 1989 reported that there is a large den- 
sity contrast between the North and South. In agreement 
with their results we find a (post shock) electron den- 
sity of (1 - 10) 10'5di7p"^V-3 for the SW Soft region and 

0.2 lO^d^^f m"3 for the NE region (d^pc is the distance 
in kpc). The range in densities for the SW Soft region 
comes from the uncertainty in the emitting volume, which 
can be anything within a factor of 10 of lO^'^dj^p^m'^. The 
estimate for the NE region was based on a filling factor of 
0.25, a coverage by the Northeastern SIS pointing of 0.2 
and a total EM of 5 10^^m~^kpc~^. The optical emission 
tells an even more interesting story: Ha observations of 
non-radiative shocks in the North give a pre-shock den- 
sity of 0.2 10^ m^'^ (Long & Blair 1990), whereas spectra 
taken from the Southwest indicate densities of more than 
Iq8 ^-3 (Leibowitz & Danziger 1983). 

The density contrast may have an accidental cause, 
but for a number of reasons we think that RCW 86 makes 
a strong case for being the result of a cavity explosion. 
In particular we think that the cavity wall has not been 
reached by the shock wave all over the remnant. With a 
cavity we mean here a tenuous region in the ISM blown 
by one or more stellar winds (Weaver et al., 1977), such 
a cavity is surrounded by a denser region, the cavity wall, 
consisting of swept-up material. The reasons for our hy- 
pothesis are: (1) the low density in the North, given the 
fact that RCW 86 is situated near the galactic plain; (2) 
the faint emission outside the main shell in Fig. |l|, most 
notably in the Northwest, implying that the wall of the 



cavity has not been reached in all directions (see Fig. ||); 
(3) Fig. 3 in Claas et al. (1989) shows that the radial X- 
ray profile of the Northeastern part of the remnant has too 
much emission coming from the shell of the SNR as com- 
pared to the Sedov model; (4) the more or less egg-shape 
reminds one of a cavity made by a star which was mov- 
ing from North to South with respect to the interstellar 
medium (Weaver et al. 1977, see their Fig. 7), or alterna- 
tively the cavity may have been blown by several stars; 
(5) the gap in the Northern shell recalls the breakouts 
seen in numerical simulations of cavity SNRs by Franco 
et al. (1991). This is even more striking in the Ha pic- 
ture by Smith (1997) where ahead of the gap there is an 
enhancement in the Ha emission; (6) a recent encounter 
of the blastwave with the cavity wall makes it more likely 
that complete electron heating has not yet occured and 
may also explain why the n^t values are rather low: most 
of the matter has only recently been shocked when the 
blastwave hit the cavity wall. 

In order to clarify point 2 we made a sketch of the 
structure of RCW 86 as we envision it (Fig. |). Note that 
the figure gives a side view and that we do not know the 
inclination, so it might also be the Southwestern region 
that is closer to us. It is clear from our sketch that we do 
not know the exact shape of the cavity, because we only 
see those parts of the cavity wall that are lit up by the 
blastwave. However, it is also clear that if the cavity has 
an egg-like shape, the remnant will also have an egg-like 
shape, but possibly less elongated. Fig. ^ also clarifies the 
relativily low emission from the center of RCW 86 (point 
3): along the line of sight the cavity wall has not yet been 
reached by the blastwave. 
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4. Conclusions 

We have presented an analysis of ASCA spectra of 
RCW 86, which wc modeled by NEI models. In general 
these models fit the data rather well, but indicate an un- 
derabundance of metals. On the other hand, we need over- 
abundances to explain the fluorescent emission of Ar and 
Fe in some parts of the remnant. The fitted models indi- 
cate a very low Uet value. 

Wc suggest that the peculiarities of the spectra may 
be the result of a non-Maxwellian electron distribution. 
If the bulk of the electrons has a very low temperature, 
most of the ions should be in a low ionization stage (e.g. 
the Li- like stage), whereas a supra-thermal tail may domi- 
nate the X-ray continuum and explain the Ar-K and Fe-K 
fluorescent emission. Note that our density calculations 
underestimate the true densities if the electron distribu- 
tion is non-thermal. 

Other SNRs may exhibit similar effects, although per- 
haps to a lesser extent. In that respect it is interesting 
that Miyata et al. (1994) and Miyata (1996) report simi- 
lar underabundances for the Cygnus Loop. 

We suggest that RCW 86 is the result of a cavity explo- 
sion. The morphology of the remnant is a good indication 
for that and oven suggests that the cavity may have been 
blown by more than one star. This is interesting in view of 
the new distance estimate by Rosado et al. (1996), placing 
RCW 86 at the distance of an OB association at 2.5 kpc 
(Westerlund 1969). If the association is real, RCW 86 is 
probably the result of a SN type II explosion. 
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